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A statistically-selected Chandra sample of 20 galaxy clusters - II. 
Gas properties and cool-core/non-cool core bimodality 
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ABSTRACT 

We investigate the thermodynamic and chemical structure of the intracluster medium (ICM) 
across a statistical sample of 20 galaxy clusters analysed with the Chandra X-ray satellite. 
In particular, we focus on the scaling properties of the gas density, metallicity and entropy 
and the comparison between clusters with and without cool cores (CCs). We find marked 
differences between the two categories except for the gas metallicity, which declines strongly 
with radius for all clusters (Z oc outside ~0.02r5oo. The scaling of gas entropy is 

non-self-similar and we find clear evidence of bimodality in the distribution of logarithmic 
slopes of the entropy profiles. With only one exception, the steeper sloped entropy profiles are 
found in CC clusters whereas the flatter slope population are all non-CC clusters. We explore 
the role of thermal conduction in stabilizing the ICM and conclude that this mechanism alone 
is sufficient to balance cooling in non-CC clusters. However, CC clusters appear to form a 
distinct population in which heating from feedback is required in addition to conduction. 
Under the assumption that non-CC clusters are thermally stabilized by conduction alone, we 
find the distribution of Spitzer conduction suppression factors, /c, to be log-normal, with a 
log (base 10) mean of -1.50 ± 0.03 (i.e. fc = 0.032) and log standai'd deviation 0.39 ± 0.02. 
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1 INTRODUCTION 

The majority of baryons in collapsed massive haloes re 
side in a hot phase, in the form of a gaseous intraclus 
ter medium (ICM), wit h the remainder predo mi nant! 
locked up in stars jpukupta. Hoean & Peebles! 1 1991 
[Gonzalez. Zaritskv & Zabludofd I2007D. This hot gas serves as 



a reservoir of material to fuel not only star formation, but also 
black hole growth, as the ultimate endpoints of radiative cooling. 
Both these processes in turn give rise to feedback from subsequent 
supernova winds (Katz 1992; Strickland_&. Stevens 2000 ) and out- 
bursts from active galactic nuclei (AGN; see lMcNamara & NulsenI 
|2003, for a recent review), respectively. 

Given the considerable effectiveness of radiative cooling in 
depleting the hot gas, its dominance of the baryon budget in cluster 
haloes emphasizes the importance o f feedback in order to restrict 
the excessi ve growth of galaxie s (e.g. lColell99 ll) and avoid a 'cool- 
ing crisis' tealogh et alj|200ll) — which has long plagued cosmo- 
logical simulations in which the effects of non-gravitational heat- 
ing are neglected (e.g. lKatz & Whit3ll993l ; ISuginohara & OstrikeJ 
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Il998l) . The same feedback mechanism(s) may also be responsible 
for ar resting gas cooling in dense cluster cores (e.g. lPeterson et al" 
'2001"), where the development of a classical 'coo ling flow' jpabianl 
1 1994) appears to be truncated dPeterson & Fabiar>i2006,) . 



A further indication of the importance of cooling is the dis- 
covery of short gas co oling times in the itmer regio ns of even 
non-cool core clusters jSanderson. Ponman & O' Sulli van 200^, 
as well as the quasi-universalit y of cooling time profiles across a 
wide range of cluster masse s l lVoigt&Fabianll2004 l Bauer et al.l 
I2OO5I ; ISanderson et alj |2006|) . despite the clear differences be- 
tween the temperature profiles of cool-core and non-co ol core 
cluste rs (e.g. Sanderson et al. 2006; Zhang et al. 2006; P ratt et al.l 
I2OO7I) . The explanation for this dichotomy in the cluster popula- 
tion is the subject of current debate (e.g. iMcCarthv et al.l 1200 8l ; 
IGuo. Oh & Ruszkowskil l2008h . but is likely to involve galaxy 
feedback, given that cosmological simulations appear to over- 
predict the abun dance of cool-cores in the cluster population (e.g. 
iKav et all |2007|) and that there is some question ove r the effec- 
tivene ss of merging in permanently erasing cool cores jPoole et"al] 
l2006h . However, the lack of evidence for strong shock heating from 
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most AGN outbursts in cluster cores ( iMcNamara & Nulsenll2007h 
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presents a challenge to understanding how feedback alone is suffi- 
cient to offset cooling losses. 

It is clear that a more complete picture of the thermodynamic 
state of hot gas in clusters is needed, in order both to solve the cool- 
ing flow problem and tackle the broader issue of feedback between 
galaxies and the intracluster medium. In pursuing these goals it is 
necessary to map the thermodynamic state of the ICM across a wide 
mass range, including both cool core and non-cool core clusters, the 
latter of which are known to be under-represented in such detailed 
studies. This enables the gas entropy to be determined, which is 
a very sensitive probe of non-gravitational processes, as well the 
metallicity, which acts as a tracer of supernova enrichment and gas 
mixing. We aim to do this using Chandra observations of nearby 
clusters, which is the only instrument able to reliably probe core 
gas properties on kpc scales, where the effects of baryon physics 
ar e greatest. Th e basis for this investigation is the statistical sample 
of ISanderson e t al. (2006, hereafter Paper I) comprising 20 galaxy 
clusters, in order to provide a representative survey of detailed clus- 
ter properties in the local universe. 

Throughout this paper we adopt the following cosmological 
parameters: Ho — 70 km s~^ Mpc~^ , Jim — 0.3 and JIa = 0.7. 
Throughout our spectral analysis we h ave used XSPEC I1.3.2t, in - 
corporating the solar abundance table of lGrevesse & Sauva 
which is different from the default abundance table. Typically this 
results in larger Fe abundances, by a factor of ~1.4. All errors are 
la, unless otherwise stated. 



2 SAMPLE SELECTION 

The objects studied in this paper comprise the statistica l sample 
of 20 galaxy clusters observed with Chandra presented in lPaper j 
The sample contains the 20 highes t flux clusters draw n from the 
63 cluster, flux-limited sample of llkebe et ^ j2002l) . excluding 
those objects with extremely large angular sizes (the Coma, For- 
nax and Centaurus clusters), which are difficult to observe with 
Chandra owing to its limited field of view. The llkebe et al.l flux- 
li mited sample was itself const ructed from the HIFLUGS sample 
of iReiprich & Bohringeil ilOOA additionally selecting only those 
clusters lying above an absolute galactic latitude of 20 degrees and 
located outside of the Magellanic Clouds and the Virgo Cluster re- 
gions. 

2.1 Re-analysis of Chandra data 

Since the original analysis of the statistical sample data in lPaper 3. 
all but one of the Chandra observations have been reprocessed 
with uniform calibration by the Chandra X-ray Center (CXC) and 
we have reanalysed all the data accordingly. The only exception is 
Abell 2256 (ObsID 1386), which was observed at a frontend tem- 
perature of - HOC: thus far, only datasets observed at -120 C have 
been reprocessed, in order to provide a uniform calibration of most 
of the data in the Chandra archiv^B Despite this, there is no indi- 
cation that our results for Abell 2256 are anomalous in any way. 

For the clusters Abell 401, Abell 496, Abell 1795, Abell 2142 
and Abell 4038, longer observations are now available and, in the 
case of Abell 1795 & Abell 2142 these are with ACIS-I, giving 
th e advant age of a wider field-of-view over the datasets analysed 
in lPaper 3 . For Abell 478, however, we retain our original ACIS-S 

^ for details, seel http://cxc.harvard.edu/cda/repro3.html| 



Name 


Obsid" 


Detector'' 


Datamode'^ 


Abell 401 


2309 


I 


F 


Abell 496 


4976 


S 


VF 


Abell 1795 


5289 


I 


VF 


Abell 2142 


5005 


I 


VF 


Abell 4038 


4992 


I 


VF 



Table 1. Clusters from the statistical sample for which new observa- 
tions have been analysed. "Chandra observation identifier. ''Denotes either 
ACIS-I or ACIS-S. '^Telemetry data mode (either Faint or Very faint). 

analysis (ObsID 1689) in favour of a newer ACIS-I observation, as 
the latter exposure is much shallower. Details of the new datasets 
analysed are given in Table[T]and key properties for the full sample 
are listed in Table |2| The dat a analysis and reduction were per- 
formed as detailed in IPaper ll using version 3.4 of the standard 
software — Chandra Interactive Analysis of Observations (CIAC0), 
incorporating CALDB version 3.4.2. 



3 DATA ANALYSIS 

Spectr al fitting was performed as described in ISanderson et al.l 
( l2006l) , using weighted response matrix files (RMFs) generated 
with the CI AO task 'MKACISRMF' for all cases, except Abell 2256 
- the only non-reprocessed dataset - where the older task 'MKRMF' 
was used instead. Spectra were fitted using an absorbed (WABS 
XSPEC component) APEC model over the energy range 0.5-7.0 keV 
for observations made with the ACIS-S detector, and 0.7-7.0 keV 
for those made with ACIS-I (as indicated in column 3 of Table|2]l. 
Spectra were grouped to a minimum of 20 counts per bin and fitted 
using the statistic. 

3.1 Cluster mean temperature and fiducial radius 

Measurements of mean temperature, T, and fiducial scaling radii 
are very important in scaling studies, in order to permit fair a com- 
parison of properties across a wide range of cluster mass. In PaperJ 
we used an iterative scheme to determine both a core-excluded 
mean temperature and rsoo (the radius enclosing a mean overden- 
sity of 500 with respect t o the critical dens i ty of the Universe), 
via the M — Tx relation of Finoeuenov et "d] (l200ll) . However, we 
have subsequently refined this method in the following two ways, 
to improve the reliability of our results. 

Firstl y, we have adopt ed the newer, Chandra-deiived M — Tx 
relation of IVikhlinin et al.l ( 2006) which is based on clusters with 
high quality observations, allowing direct measurements of gas 
properties at rsoo. Secondly, we have now excluded a larger cen- 
tral region of the data (O.ISrsoo), to remove more completely any 
contaminating emission from strong central cooling that may occur 
in the core. Our chosen annular extraction region therefore spans 
the range 0. 15-0.2r5oo; we use 0.2r5oo to exclude outer regions 
where data incompleteness begins to affect our sample, due to the 
restricted Chandra field-of-view. None the less, one advantage of 
using an outer annulus of 0. l-0.2r5oo is that it provides a bet- 
ter measure of any central temperature drop (within 0.1 rsoo), as 
this range typically brackets the peak around the cool core radius. 
Therefore we ret ain the o bjective classifications of cool-core sta- 
tus determined in lPaper j (and listed in TableO, which were based 

^ http://cxc.harvard.edu/ciao/ 
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on the significance of the temperature difference between the spec- 
trum extracted in the range 0.1-0.2r5oo compared to that measured 
inside 0.1 r 500. 

The temperatures in the M — Tx relation of IVikhlinin et al.l 

( |2006[) were based on an aperture extending out to rsoo, exclud- 
ing the central 70kpc. Since cluster temperature profiles generally 
decline monotonically with radius ou tside of any cool core (e.g. 
IVikhlinin et al.ll2005l : |Pratt et alj|2007l) . it is therefore necessary to 
allow for a bias factor that would shift temperature measurements 
within 0.2r5oo higher compared to those made within rsoo. This 
was done by c omparin g the temperatures obtained using the itera- 
tive method of 
the sample of 



Paper Bin the range 0. 15-0.2r5oo for? clusters from 
Vikhlinin et al] | |2006|) which were analysed in this 



work and in a companion analysis of galaxy groups (O' Sullivan 
et al. in prep.), namely Abell 262, 478, 1795, 2029 and MKW 4, 
MS1157 and USGCS152. The aim of this process was to deter- 
mine a multiplicative factor to apply uniformly to the temperatures 
measured in the rang e 0.15-0.2r5oo t o bring them into line with 
those measured in the I Vikhlinin et al] aperture. The determination 
of this rescale factor was also iterative, and the procedure was as 
follows. 

Firstly, a mean temperature { T) was measure d iteratively in 
the range 0. 15-0.2r5oo using the Ivikhlinin et all M — Tyi rela- 
tion to achieve convergence for each of the 6 clusters common to 
both samples. Then a rescale factor was calculated as the mean 
ratio of these tempera tures to the corresponding ones quoted in 
Ivikhlinin et alj j2006l) . Secondly, the determination of T in the 
range 0.15-0.2r5oo was repeated, using this rescale factor to ad- 
just the measured temperature before using Equation[T]to calculate 
rsoo. The process was repeated until convergence of the rescale 
factor, /, whose final value was determined to be 0.96, i.e. temper- 
atures meas ured in our aper ture were 4 per cent hotter than those 
measured by I Vikhlinin et al.L For a cluster of redshift, z, the radius 
is given by 



rsoo = 



/ X 484.7 X r"''"^ 



E{z) 



kpc 



(1) 



where / = 0.96 is a rescale factor to convert temperatures mea- 
sured in the range 0.15-0. 2r5oo with those in the range 70 kpc- 
rsoo, and 



E(z) = {l + z), 1 + {Z Qm) + 



(1 + ^) 



(2) 



A com parison between the se new temperatures and those from 
the original llkebe et al ] l l2002h parent sample is shown plotted in 
Fig-El There is sli ghtly better agreement between our new tempera- 



tures and those of 



" Ikebe et aU amongst the hott est clusters than was 
obtained using the measurements from lPaper j. alth ough the differ- 



ence betw^en_this_glot and figure 2 from Paper ll is fairly smal|f| 
The llkebe etZH2002h mean temperatures were based on a two- 
component spectral fit to the full cluster emission, which would 
likely result in lower T values for the hottest clusters, as pointed 
out in IPaper 3 . An additional factor is a possible bias in tempera- 
ture estimates for hotter clusters (>4-5 keV) resu lting from errors 
in the Chandra response matrix (as described in 'Sun et alj|2008l 
and references therein), which could lead to overestimates of T. 
Our final T values are listed in Table |2l together with the mean 
metallicity measured in the same aperture (see Section |431 >. 



^ NB the cool-core status of the points in figure 2 of jSanderson et al]j20O6l) 
was incorrectly labelled. 



-I- cool core 

O non-cool core 



O 



o 



o 
<0 



2 5 
Temperature (keV) from this work 

Figure 1. A comparison o f the new mean te mperatures from this work 
(0.15-0.2r5oo) and those of llkebe et afl J2002h . showing the line of equal- 
ity. The dots and solid lines indicate the marginal medians and interquartile 
ranges. 



Although the M - Tx relation of l Vikhlinin et alj ( l2006l) com- 
prises clusters hotter than ~2.3 keV, we none the less have verified 
that the masses of our cooler galaxy groups (O' Sullivan et al. in 
prep.) are consistent with it. Furtherm ore, we note that the recent 
Chandra study of 40 galaxy groups bv lSunetal.1 ( 12008 *) finds only 
weak evidence (at the 1.5cr level) for a steepening of the A/ — Tx 
relation in groups compared to clusters. This therefore justifies our 
inclusion of a mixture of clusters an d groups in the 7 systems used 
to calibrate our modification of the IVikhlinin et alj M — Tx rela- 
tion for temperatures measured in the aperture 0.15-0.2rsoo. 



3.2 Spectral profiles and deprojection analysis 

The results which f ollow ar e based on the deprojection analysis 
method described in lPaper i[ to derive three-dimensional gas tem- 
perature and density profiles. To stablize the fitting, the Galactic 
absorbing column and gas metallicity were fixed at values obtained 
by fitting each annulus separately prior to the deprojection. Con- 
sequently, the gas metallicity results presented in Section 1431 are 
projected quantities. For some clusters, it was necessary to freeze 
the absorbing column at the galactic HI value, since unfeasibly low 
values were otherwise obtaine d in many of the annular bins: full 
details can be found in lPaper j . Similarly, in a number of cases the 
deprojected temperature had to be fixed a t its pro jected value, to 
produce a stable fit, exactly as was done in lPaper j . 



4 RESULTS 
4.1 Gas density 

Fig. |2| shows the gas density as a function of scaled radius for the 
sample, colour-coded by mean temperature and s plit according to 
cool-core status. Following lSanderson et al.l ( l2006l) . in order to clar- 
ify the underlying trend in each profile, the raw data points have 
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Name ObsID" Detector'' Redshift Mean kT'^ 'fsoo Mean metallicity'^ Cool-core status Si -^^ Index'' 











(keV) 


(kpc) 


(Solar) 




(keV cm^) 




NGC 5044 


798 


S 


0.008 


1.17+°°! 

— 0.05 


512+1° 

— 10 


0.37+°-°° 
—0.06 


cc 


50+?" 

— 17 


0.71+°-3° 
—0.09 


Abell 262 


2215 


S 


0.016 


— 0.09 


692+1? 

— 15 


0.38+°-i2 

—0.07 


cc 


126+*» 

— 20 


0.91+°?° 

—0.13 


Abell 1060 


2220 


I 


0.012 


2.92+°-" 

— 0.11 


829+}^ 


0.50+°°^ 

—0.07 


non-CC 


191+° 

— 7 


0.44+°-°? 
—0.07 


Abell 4038 


4992 


I 


0.030 


3.04+°-°^ 

— 0.09 


840+1° 

— 13 


0.60+°°^ 

—0.07 


non-CC 


136+* 


0.55+°-?! 

—0.14 


Abell 1367 


514 


s 


0.022 


3.22+°-?8 
— 0.18 


869+?5 

""'^ — 26 


0.34+°i° 

—0.09 


non-CC 


275+" 

— 31 


0.28+°-i? 
—0.31 


Abell 2147 


3211 


I 


0.035 


3.69+°" 

— 0.18 


— 23 


0.28+°ll, 

—0.10 


non-CC 


281+^^ 

— 59 


0.43+°-?I 
—0.16 


2A 0335+096 


919 


s 


0.035 


4.09+°" 

— 0.13 


980+1? 

^"" — 16 


0.79+°ll 

—0.11 


CC 


99+?° 

— 20 


1.12+°-?2 

—0.16 


Abell 2199 


497 


s 


0.030 


4.50+°-?° 

— 0.24 


I033+?a 

v-^^ _2g 


0.79+°ll 

—0.14 


CC 


165+00 

— 33 


0.76+°°° 

—0.06 


Abell 496 


4976 


s 


0.033 


4.80+°-" 

— 0.14 


1067+18 

— 16 


0.66+°°^ 

—0.08 


CC 


162+^5 

— 34 


0.93+^'^^ 

—0.07 


Abell 1795 


5289 


I 


0.062 


5.62+°-36 

— 0.35 


1144+^1 

— 38 


0.26+°l° 

—0.10 


cc 


— 35 


0.99+n n2 
—0.08 


Abell 3571 


4203 


s 


0.039 


6.41+°?? 

— 0.23 


1239+0? 

— 23 


0.75+°-ll 

—0.11 


non-CC 


279+° 

— 7 


0.44+°-°l 

—0.04 


Abell 2256 


1386 


I 


0.058 


6.52+0-39 

— 0.36 


1239+?^ 

— 36 


0.98+°-?? 

—0.21 


non-CC 


344+?t 
— 158 


0^47+°-'?I 

—0.16 


Abell 85 


904 


I 


0.059 


6 64+°-2° 


125lt« 


0.561°:°? 


CC 


1Q3+I6 


90+0.14 

"■^"-0.09 


Abell 3558 


1646 


s 


0.048 


7 ,7+0.49 
'■"-0.46 


1309t« 


1.001°:- 


non-CC 


304+^3 


62+°-3° 

"■"^-0.11 


Abell 3667 


889 


I 


0.056 


7 60+°-^* 

'■""-0.37 


13451^4 


5J+0.10 

"■-"-0.10 


non-CC 


407lf, 


57+0.17 
".-"-0.09 


Abell 478 


1669 


s 


0.088 


n 7^+0.26 

"■^-'-0.26 


1381123 


50+0.07 
" -'"-0.07 


CC 


190+1* 


1 03+°-2l 
i."-'-0.06 


Abell 3266 


899 


I 


0.055 


o ,n+0.67 

».3»_o.43 


1/17+58 
''••"-38 


39+0.11 

"■^^-0.11 


non-CC 


528119 


f, 40+0.09 

0.48_Q Qg 


Abell 2029 


4977 


s 


0.077 


8 96+°-3° 


1452+25 


"■60-0.07 


CC 


9<:7 + 20 
^J'_49 


90+0.21 

"■^"-0.08 


Abell 401 


2309 


I 


0.074 


^■'°-1.06 


1471111* 


26+0.22 
"■^"-0.23 


non-CC 


477+55 


45+0.15 
"■^^-0.11 


Abell 2142 


5005 


I 


0.091 


9 50+° '*3 

^•^"-0.42 


14871^^ 


44+0.07 

"■^^-0.07 


non-CC 


2951^3 


94+0.06 
"■^^-0.17 



Table 2. Key properties of the sample, listed in order of increasing temperature. "Chandra observation identifier. ''ACIS detector. '^Measured between 0.15 
and 0.2r5oo (see Section lxTt . '^Parameters of the power-law fit to the entropy profile (see Section l4.2.H . Errors are la. 



been fitted with a locally weighted regression in log-log space, us- 
ing the task 'LOESS' i n version 2.7 of the R statistical software 
environment packag^B i IR Development Core Teaml2008h . It is im- 
mediately apparent from Fig. [2] that the profiles do not scale self- 
similarly, which would imply constant density at a given fraction 
of rsoQ. There is a wide dispersion in the core, with a factor of 
~30 range in pgas at O.Olrsoo, dropping to a factor of ~5 spread 
at 0. IrsoQ. At larger radius ( beyond the peak in the temperature 
profile at ^0. ISr.g^no. IPaper^ there is a clear convergence in the 
profiles, despite the diminishing data coverage due to the limits of 
the Chandra field of view. 

It is clear from Fig. [2] that the cool-core (hereafter CC) clus- 
ters have systematically denser and more cuspy cores, c onsistent 
with the decline in gas temperature (scaling with r^O'*; [p^erj) 
that must be counteracted by a rising density in order to main- 
tain pressure equilibrium. The near power-law shape of the density 
profiles in the CC clusters is also consistent wi th a cooling domi- 
nated regime, as indicated by the simulations of lEttori & Brighentil 
( |2008[) . who find pgas oc r~^'^ for evolved cool-core clusters (after 
10 Gyr). Within each cool-core category, there is also a systematic 
trend towards lower ICM densities in cooler systems. This may re- 
flect gas depletion due to cooling out of the hot phase, or could 
be caused by non-gravitational heating expelling material out of 
the core. The se findings are consistent with the recent study of 
ICroston et al. II2OO8), based on a representative sample of 3 1 clus- 
ters analysed with XMM-Newton. 

4.2 Gas entropy 

The entropy of the gas is a key parameter, which provides a mea- 
sure of the thermodynamic state of the ICM and is conserved in 

* http://www.r-project.org 




2 4 6 8 
Temperature (keV) 



0.01 0.1 
Radius (rsoo) 

Figure 2. Gas density profiles for each cluster, scaled to rsoo and coloured 
according to the mean cluster temperature, depicted by the inset histogram. 
Each curve represents a locally weighted fit to the data points, to suppress 
small-scale fluctuations (see the text for details). 



any adiabatic process (see iBowej 1 19971; iTozzi & Noiiiiai] I2OOII; 
IVoit et al.ll2002l for example). We define entropy as 5* = kT/pgi^ 
(e.g. I Ponman. Cannon & Navarrclll99 9'). which implies S cc kT 
for self-similar clusters. Within a virialized halo, the gas entropy is 
initially set by shock heating, which leads to a r adial variation of 
the fo rm 5* oc ri 1 for a simple spherical collapse jTozzi & NormanI 
l200lh . Despite radiative cooling in the cores of clusters lower- 
ing the entropy, its power-law variation with radius is nevertheless 
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1000 



100 



10 



cool core 
non-cool core 




Temperature (keV) 



10 100 
Radius (kpc) 

Figure 3. Gas entropy profiles as a function of physical radius, colour- 
coded by cluster mean temperature. Each curve represents a locally 
weighted fit to the data points, to suppress small-scale fluctuations (see the 
text for details). 



largely preserved, wit h approximately the same logarithmic slope 
as in the outskirts (e.g. lEttori & Brighentill2008h . leading to the ex- 
pectation of a near-proportionality between entropy and physical 
(i.e. unsealed) radius, regardless of halo mass. This can be seen in 
Fig. [3] which shows gas entropy profiles as a function of radius in 
kpc for the sample, separated by cool-core status and colour-coded 
by mean temperature. 

The distinction between CC and non-CC clusters is very clear 
from Fig. |3] with the latter having significantly higher entropy in 
the core and the divergence between the two types occurring within 
~40kpc, at an entropy level of ^SOkeVcm"^. The CC entropy 
profiles are tightly grouped and show no obvious sign of the tran- 
sition between the shock heating regime in the outskirts and the 
cooling-dominated core. In contrast, the non-CC profiles show a 
much larger scatter in normalization, which may reflect the diver- 
sity of heating processes affecting them. Alternatively, it could be 
that cooling acts to regulate the entropy of CC clusters so as to 
suppress the scatter between them, which may also account for the 
apparent universality seen in the cooling time profiles of the sample 
jPaper ih . 



4.2.1 Entropy profile fitting 

In order to characterize the form of the entropy profiles plotted in 
Fig. [3] we fitted power-laws to each cluster profile to quantify the 
logarithmic slope of the relationship. However, it can be seen from 
Fig.[3]that at small radii the profiles tend to flatten and therefor e de- 
viate from a simple power law form (see also Donahue et al. 2 0061 . 
for example). Nevertheless, a power law provides a simple and rea- 
sonably effective description of the majority of the profile in all 
cases which enables the 'flatness' of the curves to be character- 
ized. We have additionally employed a quantile regression tech- 
nique to perform the fitting, in order to provide resistance to any 
such outliers. This form of regression minimizes the sum of abso- 
lute residuals, rather than the sum of squared residuals, and thus is 
analogous to de termining the median as an estimator of the mean 
( lKoenker|[2003) . The algorithm used is the 'RQ' function from the 
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Figure 4. Residuals from the power law fits to the entropy profiles, nor- 
malized by the predicted value and plotted against scaled radius for each 
cluster. A single value is omitted (0.76 at 0.0047r5oo, for the innermost 
point of 2A 0335-1-096) to optimize the y-axis scale range for viewing the 
data. 



QUANTREG package in p[f| and was executed as an unweighted lin- 
ear fit in log-log space. Table |2] lists the best-fit normalization (at 
0. Irsoo) and power-law index (i.e. logarithmic slope) for each clus- 
ter, together with the corresponding la errors. 

To assess the suitability of a power law for describing the en- 
tropy profiles, we show in Fig. |4] the residuals from the best-fit 
model as a function of scaled radius, for each cluster. We have nor- 
malized the residuals by the predicted model values rather than the 
measurement errors on the entropy, since the fit was performed us- 
ing unweighted quantile regression rather than by minimizing x^- 
This approach allows for the fact that real clusters exhibit intrin- 
sic deviations from simple radial models in excess of the statistical 
scatter associated with measurement errors. It can be seen that in 
general the power law fits do a reasonable job of describing the 
data, with the majority of residuals contained within ^5 per cent of 
the best fit. However, while there is no obvious sign of any strong 
systematic trends with radius, there is some indication of an ex- 
cess entropy above the model at very small radii (<0.03r5oo) for 
a few of the clusters, consistent with the fl attening in the entrop y 
profiles of nine cool core clusters found bv lDonahue et al.l ( I2OO6I) . 
Nevertheless, we restrict our use of the power law fits to the en- 
tropy profile to a larger radius (0. Irsoo), which is not affected by 
any such systematic departure from a power law form. 

There is no evidence of any systematic variation in the loga- 
rithmic slope of the entropy profile with mean temperature across 
the sample, as plotted in Fig. |5] However CC and non-CC clus- 
ters are clearly segregated in the plot, occupying higher and lower 
values of the entropy index, respectively. The only exce ption is 
Abell 2142 - the archetypal 'cold front' cluster (Markevi tch et al.J 
I2OOOI) - which was classified as a non-CC cluster in lPaper A as its 
cool co re ratio was found to be only marginally significant (~ 2cr; 
IPaper ih . Notwithstanding the cool core status of the clusters, it is 
reasonable to ask if there is any evidence for bimodality in the dis- 
tribution of logarithmic slope values in Fig. [5] and we address this 
specific issue in Section[5] 



See the tutorial at |http://www.astrostatistics.psu.edu/su07/R/reg.html| 
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Figure 5. Best fit logarithmic slope of the entropy profile as a function 
of mean system temperature. The dashed line marks the cross-over point 
between the two Gaussian distributions of values from Section \5\ and the 
outlier Abell 2142 is labelled (see text for details). 



4.2.2 Entropy scaling 

The effectiveness of entropy as a probe of non-gravitational pliysics 
in tlie ICM is best exploited by studying its variation with mean 
temperature, which is expected to be hnear in the case of simple 
self-similarity. However, the strong dependence of entropy on ra- 
dius, raises the issue of at what point to measure entropy in order 
to consider its variation from cluster to cluster. The issue is further 
complicated by the fact that the gas entropy within a cluster bridges 
two distinct regimes: the cooling-dominated core and the outskirts, 
where the effects of shock-heating from infall prevail. The first 
measureme nts of entropy scaling were made at a fiducial radius 
of 0. 1r2oo jPonman et alj[l999l : iLlovd-Davies. Ponman & CanonI 
120001) in order to sample the region between these two regimes. 
However. IPratt. Amaud & Pointecouteaul jlOOfil) find consistent re- 
sults for the scaling with temperature of entropy when measured 
a range of radii at fractions of r2oo of 0.1, 0.2, 0.3 and 0.5 in the 
range 6 = 0.5-0.6 for S oc t'. 

Taking the power-law fits to the entropy profiles performed 
above we plot in Fig. |6] entropy vs. T for the sample using the 
normalization value, corresponding to a fiducial radius of O.lrsoo- 
Also shown are power-law fits to the CC and non-CC clusters 
separately, performed in log-log space using the BCES weighted 
orthogonal regression method of Akritas & Bershady ( 1996). The 
logarithmic slopes of these lines correspond to an entropy scaling 
of the form S (x t'' with 6 = 0.66 ± 0.10 and 0.71 ± 0.21 for 
the CC and non-CC clusters, respectively. These results are in good 
agreement with the value of 0.65 ± 0.05 similarly determined by 
IPonman. Sanderson & Finoguenovl ( l2003h . at a fiducial radiu s of 
0.1r2oo, as well as with the values found bv lPratt et alj ( |2006|) at a 
series of radii and using a number of regression methods. 

The CC and non-CC points are generally well separated with 
respect to their corresponding best-fit power laws in Fig. |6l How- 
ever, in the case of the CC cluster Abell 262, the (albeit large) error 
bar overlaps significantly with the non-CC best-fit line. While this 
poor cluster does have a strong negative central temperature gradi- 
ent, it is also known to possess a number of cavity and ripple fea- 
tures coincident with low frequency radio emission jBlanton et al.l 
l2004h . It is therefore possible that AGN activity in the core might 
have boosted the entropy at O.lrsoo to shift this cluster towards the 
non-CC regression line. Nevertheless, a number of the other CC 
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Figure 6. Entropy at O.lrsoo as a function of mean temperature, from a 
power-law fit to the entropy profile of each cluster. The lines are the best- 
fit power-law, with logarithmic slopes of 0.66 ± 0.10, 0.71 ± 0.21 and 
0.92 ± 0.12 for the CC, non-CC and combined clusters, respectively. 



clusters also host ghost cavities (see table 2 in lPaper j) , which evi- 
dently have not significantly biased their location in Fig.|6] 

While the separate fits to the CC and non-CC clusters in Fig.|6] 
yields results consistent with a no«-self similar entropy scaling, the 
two sets of points are clearly offset from each other. It is interesting 
to note that a BCES orthogonal regression fit to all the 20 clusters 
combined yields a much steeper logarithmic slope of 0.92 ± 0.12, 
which is consistent with self-similar scaling. To check this result, 
we repeated the individual power-law fits using a pivot-point of 
O.lSrsoo, and performed the same regression of the normaliza- 
tion versus mean temperature to obtain the following logarithmic 
slopes: 

1.06 ± 0.16 (all clusters combined) 
0.65 ±0.10 (CC) 
0.65 ± 0.26 (non-CC). 

Thus it is possible that similarity breaking can occur in the sepa- 
rate CC/non-CC cluster sub-populations in such a way as to pro- 
duce fully self-similar scaling in an analysis which combines the 
two types. On the other hand, if we exclude the coolest system 
(NGC 5044) from the fit, the resulting best-fit slope for the whole 
sample (0.89 ± 0.22) is consistent at the la level with the values 
obtained for the separate CC and non-CC sub-samples. 

Given the broad agreement in the dependence on temperature 
of the entropy scaling, we show in Fig.|7]entropy profiles a s a func- 
tiqii of r / rsoo normalized by the factor y" *'^ from Ponm an et al.l 
120031) . It can be seen that this empirical scaling brings the curves 
into fairly close alignment. There is, however, some indication of 
a systematic dispersion in the profiles, indicating that gas entropy 
may vary less strongly with cluster mean temperature than yO-^^^ 
also suggested by the results of Pratt et al. ( 2006) when measuring 
entropy at 0.1r2oo- In order to explore the entropy scaling across a 
range of radii, rather than at at particular spot value, we modelled 
the variation of entropy in terms of both scaled radius and mean 
temperature using the following expression: 
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Figure 7. Entropy profiles, scaled with the lPonman et alj )20Q3|) empirical 
factor of T"'^^ vs. scaled radius, colour-coded by cluster mean temperature. 
Each curve represents a locally weighted fit to the data points, to suppress 
small-scale fluctuations (see the text for details). 



Sample 


S" ( keV cm^) 


a 


h 


All 
CC 
NCC 


259+^ 


Ogj + O.lO 

"■^-"-o.os 
1 05+°°^ 

5J+0.09 

"■-'^-o.oe 


56+"" 
"■-'"-o.os 
49+0-07 

52+0-08 

"-•^^-0.07 



Table 3. Results from fitting Equation |3] to the entire sample, as well as 
separately to the combined cool-core and non-cool core clusters. Errors are 
Icr; see text for details. 



This parametrization has the advantage of allowing the scaling of 
entropy with both radius and temperature to be determined simulta- 
neously. The results of performing an unweighted quantile regres- 
sion fit in log-log space of Equation[3]to the entire sample are sum- 
marised in Table[3] Only data in the range 0.02 < r/rsoo < 0.2 
were included, to ensure completeness in the radial coverage. 

The resulting values of a agree well with the results from the 
individual power law fits to each cluster described previously, in 
Section [4.2. II It can also be seen that the fit to the entire sample 
yields a dependency on mean temperature such that 5 cx T , with 
6 — 0.56+Q slightly lower but still in agreement with the scaling 

which was evaluated 

at a r adius of 0.1r2oo — roughly O.lSrsoo dSanderson & Ponniiai] 
l2003h . A fit to the combined CC and non-CC clusters separately 
yields values of 6 = 0.49^o;?3 and h = 0.52^°;°*, respectively, 
demonstrating that this weaker dependence of entropy on system 
temperature persists. This result is strongly inconsistent with a self- 
similar scaling of 5 « T (i.e. h — \) demonstrating the im- 
pact of non-gravita tional physics in influencing the hot gas (see 
IPotmian et al.ll20()3l for ex ample, for furt her discussion), but ex- 
tends beyond the finding of IPonman et al . in that it applies across 
a range of radii spanning an order of magnitude, as opposed to be- 
ing determined at only a fixed spot value. Furthermore, the fact 
that this result is essentially unchanged when fitting the CC and 
non-CC clusters separately (see Table [3}, demonstrates the uni- 
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Figure 8. Gas metalHcity ( Greves se & Sauva ] |l998l abundances) versus 
mean temperature, both measured in the range 0.15 ^ f/fsoo ^ 0-2. 

versality of this modified entropy scaling, at least in the range 
0.02 < r/rsoo < 0.2. 

However, the difference between CC and non-CC clusters is 
very clear when considering the variation of entropy with radius, 
5* oc r". Here the corresponding logarithmic slopes are a = 
1.05lEJ;°5 and O.SllEJ g^, respectively, compared to 0.95lg;J^ for 
the entire sample. The flatness of the best-fit non-CC entropy pro- 
file compared to that for the CC clusters is further underscored by 
the difference in normalization at the fiducial 'pivot point' in Equa- 
tion[3](at O.lrgoo for T = 5 keV), yielding values of 5" = 259+4 
keV cm^ compared to IGII^q keV cm^, respectively. The results 
from the combined fit of Equation [3] are in good agreement with 
the results from the power-law fits to the individual cluster entropy 
profiles. T hey are also consiste nt with the log arithmic slope of 
0.95±0.02 jPiffaretti et alj2005h and 1.08±0.04 jPratt et alj20o3) 
from two separate XMM-Newton analyses of cool-core clusters, as 
well as the values of slopes in the range 1-1.3 fr om the Chandra 
analysis of 9 CC clusters bv lDonahue et al.l ( l2006h . 

4.3 Metallicity of the ICM 

The metallicity of the ICM is an important tracer of galaxy feed- 
back via supernova-driv en winds, which eject both m etals and en- 
ergy into the hot gas (e.g. lStrickland & Stevensl200(]|) . Fig.[8]shows 
the mean metallicity against mean temperature for the sample (see 
TableO, both evaluated in the range 0.15 ^ r/rsoo 0.2, iden- 
tified according to cool-core status. There is no indication of any 
trend in the data or of any systematic difference between CC and 
non-CC clusters; the mean value across the entire sample is 0.55 
Solar, with a standard deviation of 0.22. 

To explore radial trends we use the metallicity values deter- 
mined in the projected annular spectral fitting, which were taken as 
fixed inputs in the deprojection analysis as described in Section |T2l 
and Paper 1. The fact that these values are no t depro j ected will tend 
to smo oth out gradients slightly. However. iRasmussen & PonmanI 
( l2007h demonstrate that this is a small effect and that deprojected 
metallicity data typically suffer from large instabilities which add 
significantly to the noise, particularly given the poorer constraints 
on this parameter compared to temperature or density. 

Fig. [9] shows the individual projected metallicity profiles as 
a function of scaled radius for the sample, split by cool-core sta- 
tus and colour-coded by mean temperature. The curves represent a 
locally-weighted fit to the data, using the LOESS function described 
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Figur e 9. Projected gas metalUicity profiles (using iGrevesse & Sauva] 
119981 abundances) for each cluster, scaled to rsoo and coloured according to 
the mean cluster temperature, depicted by the inset histogram. Each curve 
represents a locally weighted fit to the data points, to suppress small-scale 
fluctuations (see the text for details). 

in Section 143] and it can clearly be seen that the metallicity de- 
clines with radius beyond ~0.01-0.02r5oo in almost every case, 
for both CC and non-CC clusters. Within this radius the dispersion 
in metallicity increases noticeably; some clusters show a strongly 
peaked metallicity profile (in particular Abell 85 and Abell 2029), 
whereas others have metal-deficient central cores. More specifi- 
cally, the three most peaked profiles are all hot clusters, whereas 
the coolest clusters show sharp central declines. It is not clear what 
is responsible for this strong divergence in behaviour, but it occurs 
on a scale of roughly 0.02r5oo (^20-30 kpc), which is comparable 
to the size of the central galaxy. This also coincides with the region 
where the gas density and entropy profiles also become widely dis- 
persed (Figs.|2]and|7l(. 

Regardless of any difference between CC and non-CC clus- 
ters at small radii, it is clear that the point where the metallic- 
ity begins to decline strongly towards the cluster outskirt s lies 
well within the typical cooling radius of ~0.15r5oo jPaper j) . The 
fact that the metallicity profiles of all the clusters are quite sim- 
ilar beyond this point suggests that the enrichment of the ICM 
is insensitive to those factors responsible for determining cool- 
core status. The finding that metallicity declines strongly with ra- 
dius in all the cluster s is ap parently at odds with the study of 
iDe Grandi & Molendil ( 1200 ll) , who concluded that the abundance 
profiles for the 8 non-CC clusters in their BeppoSAX sample were 
consi stent with being const a nt. Ho wever, examination of figure 2 
from lDe Grandi & Molendil ( 1200 ll) shows that, with the exception 
of only a few outliers, there is an indication of a general decline 
in metallicity with inc reasing radius in th eir non-CC clusters. Fur- 
thermore, we note that lSaldi et al.l j2007h also found no difference 
between CC and non-CC gas metallicity profiles, outside ~0. lri8o 
in their analysis of 12 hot clusters observed with Chandra. 

To highlight the trend in metallicity with radius we show the 
mean CC and non-CC profiles in Fig. (TO) grouped to a total of 5 
bins in each case, with error bars indicating the standard deviation 
within each bin. Apart from a slightly larger dispersion in the non- 
CC bins, there is very little difference between the two profiles and 
no indication in either case of any flattening off at large radii. In- 
teresti ngly, the BeppoSAX average CC profile of lOe Grandi et all 
( |2004|) shows a sharp levellin g off in metallicity at around 0.4 So- 
lar dGrevesse & Sauvall 19981 , abundances) at 0.2r2oo, which corre- 
sponds to roughly O.Srsoo - close to the limit of our Chandra data 
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Figure 10. Average gas metallicity profiles for cool-core and non-cool 
core clusters. Each bin represents the mean of a roughly equal number 
of points and the en'or bars show the standard deviation in both direc- 
tions. The dashed line is the best-fit power-law to the unbinned data beyond 
0.014r/r5oo and the dotted line is a fit to the iron abundances in 15 galaxy 
groups from Rasmussen & Ponman (2007, see text for details). The points 
and horizontal lines indicate the medians and interquartile ranges of the CC 
and non-CC raw data points. 

but otherwise consistent with the outer bin of Fig.fTO] so we cannot 

rule out a flattening beyond t his point. 

However, in contrast, the lOe Grandi et aL Uiooi non-CC pro- 
file is both flatter and lower in normalization (albeit only at ~2(t) 
and thus rather different from our own. Of their 10 non-CC clusters, 
only three are present in our sample (A 1367, A2256 and A3266; 
they classifiy A2142 and A357I as CC clusters), and their profiles 
for these systems appear to be reasonably consistent with our own, 
within the region covered by C handra. The difference s between 
CC and non-CC clusters seen by iDe Grandi etal.l ( l2004l) appear to 
originate in the somewhat low and fairly flat abundance profiles 
of several non-CC clusters in their sample, whereas the non-CC 
clusters in our sample don't appear to include s uch cases. Notwith- 
standing this difference in samples, we note that lBaldi et alj ( l2007h 
measure a similar radial decline to us, with essentially no difference 
between CC and non-CC clusters, albeit with the exception of their 
central bin. As mentioned above, cluster metallicity profiles appear 
to exhibit greater dispersion on the scale of the central galaxy, so 
some difference from sample to sample is to be expected in this 
region. 

The data in Fig. [TO] appear to be reasonably consistent with 
a power law, so we have fitted such a function in log-log space, 
using the quantile regression method outlined in Section |4.2.1| In 
order to exclude the core region where the profile flattens and the 
profile diverge substantially (r ~ 0.014r5oo; Fig.|9j, we performed 
separate fits inside and outside this radius, combining both the CC 
and non-CC data. The results for the inner region (r < 0.014r5oo) 
are 

logio(^/^Q) = 0.3912:^2 log,o(r/r5oo) + 0.821° °', 
while those for the outer region (r > 0.014r5oo) are 
logio(Z/ZQ) = -0.3llH5logio(r/r5oo) - 0.5X1^?, 
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Figure 11. Projected gas metallicity vs. entropy, with locally weighted fits 
to the CC and non-CC clusters plotted as solid and dashed lines, respec- 
tively. 



By comparison, the metallicity profiles of galaxy groups appear 
to be somewhat steeper in their decline at large radii. The recent 
analys is of 15 group s with Chandra data by Rasmussen & PonmanI 
( l2007h - also using Grevesse & Sauvall il9m abundances- found 
logarithmic slopes of -0.66 ± 0.05 and -0.44 for the combined ra- 
dial profiles of iron and silicon, respectively, across their sample. 



compared to our slope of -0.3112:21 

for the mean me tallicity (see 

Fig. not . At a radius of r.»inn ■ iRasmussen & PonmanI estimate an 
average iron abundance in groups to be ~0. 1 , whereas the extrapo- 
lation of the above best fit implies a mean metallicity of Z ~ 0.15 
in clusters at the same radius. 

The similar abundance patterns between CC and non-CC clus- 
ters can be seen by examining the gas metallicity as a function 
of entropy, plotted in Fig. [TT] Also plotted are locally-weighted 
curves for the combined points in each category, which demon- 
strate that there is very little difference between them, aside from 
the greater radial coverage in the centres of CC clusters. Since 
the gas entropy scales roughly linearly with radius (Fig. O, the 
trend between metallicity and entropy mimics that seen in Fig. |9] 
However, a striking aspect of Fig. [TT] is the almost complete ab- 
sence of low-metallicity gas (i.e. < 0.4-0.5 Solar) with low en- 
tropy (5* < 200 keV cm^). In general it appears that the most 
enriched gas has low entropy, although the smoothed local re- 
gression suggests a turnover in this trend towards lower metal- 
licity below ~30 keVcm^. Interestingly, this entropy level cor- 
responds to the threshold below which star for mation appears to 
take place in the central g alaxies of CC clusters jVoit et al]|2008l : 
iRaffertv. McNamara & Nulse n 2008). Thus, a reversal of the in- 
verse trend between entropy and metallicity may reflect the loss 
of the most enriched gas from the hot phase in fuelling such star 
formation. 



5 CLUSTER BIMODALITY AND THERMAL 
CONDUCTION 

Clusters can be divided according to the presence or absence of a 
cool core, and it is clear that the properties of these two categories 
differ substantially in terms of their temperature, density and en- 
tropy profiles. However, an important question to ask is, are these 
two types merely separate parts of the same continuous distribution, 
or do they really constitute distinct populations? If the latter holds 
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Figure 12. Kernel-smoothed (with a Gaussian of standard deviation 0.07) 
probability density plot of the best-fit logarithmic slope of the entropy pro- 
file (soHd line), showing the raw values as randomly 'jittered' points over 
the X-axis. The dashed curves are the components of the best-fit bimodal 
Gaussian distribution to the data (see text for details). With the exception 
of Abell 2142, all the clusters above 0.7 are CC and all those below are 
non-CC. 



true, then this would have interesting implications for models of 
self- similarity breaking via feedback and/or other non-gravitational 
processes. 

To address this issue, we return to the power-law fits to the 
entropy profiles described in Section [4.2. 11 The probability distri- 
bution of logarithmic slopes is plotted as a kernel density estimate 
in Fig. [12] with the positions of the raw values also indicated. It 
can be seen that two peaks are visible in the distribution, sugges- 
tive of bimodality in the cluster population. We performed a maxi- 
mum likelihood fit to the unbinned entropy slope values with both 
a single and a pair of (equally-weighted) Gaussians, using the FIT- 
DISTR function in the MASS package in R. The single Gaussian 
best-fit mean value was jj, — 0.70 ± 0.05 with standard deviation 
of (7 = 0.25 ± 0.04, and the best fit double-Gaussian values were: 



Hl = 0.48 ± 0.04 
fi2 = 0.92 ± 0.04 



crl = 0.10 ±0.03 
ct2 = 0.11 ±0.03. 



The two separate Gaussians are plotted as dashed curves in Fig.lI2l 
and align closely with the peaks in the smoothed density distribu- 
tion; the cross-over point between these two Gaussians is 0.69 and 
is plotted as the dashed horizontal line in Fig. [5] 

In order to provide a quantitative assessment of the putative 
bimodality in the pop ulation, we cal culated both the Akaike infor- 
mation criterio n (AIC, AkaikellT974h and the Bayesian information 
criterion (BIC, ISchwarzl I978|). These ar e standard statistics used 
in model selection (e.g. see Eiddlell200l for a discussion of their 
relative merits), based on the negative log-likelihood penalised ac- 
cording to the number of free parameters, with the model having 
the lowest value always being preferred. The corresponding values 
obtained are AIC = 4.65 (-0.49) and BIC = 6.64 (3.49) for the sin- 
gle (double) Gaussian model, implying a change in BIC of 4.0 in 
favour of the two component fit. Differences in BIC of between 2 
and 6 indicate positive evidence against the model with the greater 
BIC value, with values above 6 indicating strong evidence agains t 
the model with the greater BIC value (e.g. lMukheriee et al.lll998h . 
demonstrating that a bimodal distribution is clearly favoured over 
a unimodal distribution. Given that the sample was statistically se- 
lected, this therefore implies that two distinct categories of cluster 
exist. 

The two separate distributions of entropy profile logarithmic 
slopes match the cool-core classification (with the exception of 
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Abell 2142; see Section l4.2. It . show similar dispersion (a ~ 0.1) 
and are well-separated (the means differ by ~4 standard devia- 
tions). The means of the distributions are also in good agreement 
with the results obtained by fitting Equation [3] to the CC and non- 
CC clusters separately, as summarised in table[3] 

While a power law fit describes the entropy profiles of CC 
clusters quite well, it can be seen from Fig.|7]that the non-CC pro- 
files appear to flatten increasingly at small radii. This raises the 
possibility that estimating the gradient from a power law fit (al- 
beit robustly) could bias the results. Therefore, to check our con- 
clusions regarding bimodality, we have also evaluated the entropy 
profile logarithmic slope at O.OSrsoo as estimated from a locally- 
weighted fit to the data (in log-log space), using the LOESS task 
in R (c.f. Fig. [7}, with a heavier smoothing (using a value of the 
'span' parameter of 2), to increase the stability of the gradient mea- 
surements. Repeating the above analysis leads to best fit bimodal 
populations, given by 



/il' = 0.26 ± 0.09 crl' 
/i2' = 0.97 ± 0.03 cr2' 



0.28 ±0.07 
0.08 ±0.02. 



The corresponding values obtained are AIC = 25.4 (18.4) and 
BIC = 27.4 (22.4) for the single (double) Gaussian model, imply- 
ing a change in BIC of 5.0 in favour of the two component fit. It 
can be seen that the non-CC slope is indeed flatter at this smaller 
radius (although with much greater dispersion), but that the conclu- 
sion that the population is bimodal is verified. 



5.1 Conduction and thermal balance 

One of the most intriguing aspects of the ICM is its resistance to 
runaway radiative cooling. Although galaxy feedback is a plausible 
mechanism for maintaining (or nearly maintaining) thermal bal- 
ance in the centres of cool-core clusters, it is less clear if it can sim- 
ilarly affect non-CC clusters, which none th e less also have short 
cooling times in their inner regions jPaper j) . This is because non- 
CC clusters show no sign of mass drop-out necessary to fuel AGN 
outbursts or supernova winds. However, heat transport by thermal 
conduction is likely to play an important role in such cases and is 
also very effective at stabilizing cooling within CC regions, pro- 
vided the gas de nsity is not too high (e.g. lConrov & Ostrikej20o3 : 
IGuo et al.ll2008h . 

Circumstantial evidence for the effectiveness of thermal con- 
duction in counteracting cooling can be seen in the tempe rature 
profiles of CC clusters. As pointed out bv lVoigt et alj i2002l) . cool 
cores at the limit of stablization by conduction would have temper- 
ature gradients of the form T oc r" '* in the case of bremsstrahlung 
emission, flattening to T oc r"'^ for line-dominat ed emis sion, 
which concurs well with observations of CC clusters jPaper J) and 
groups ( O' Sullivan et al. in pre p.), respectively . Furth ermore, recent 
work bvl Voitetal.1 j2008l) and lRaffertv et al.l ( |2008|) indicates that 
star formation in the cores of clusters, resulting from unchecked 
cooling, only occurs when the gas entropy drops below a certain 
level (< 30 keV cm^) which matches the threshold below which 
gas becomes thermally unstable against conduction. 

The critical scale for thermal conduction is given by the Field 
length, which must exceed the characteristic size of the system in 
order to allow thermal bal ance to be m aintained by smoothing out 
tempe ratur e fluctuati ons l lFieldll965l) . Following iDonahue et al.l 
( l2005l) and lVoit et al] j2008l) . the Field length, Af, can be approxi- 
mated as 
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Figure 13. Implied Spitzer conductivity suppression factor as a function 
of scaled radius, grouped by cool-core type and plotted with identical axes. 
Each curve represents a locally weighted tit to the data points, to suppress 
small-scale fluctuations (see the text for details). 
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where k is the Spitzer conduction coefficient with suppression fac- 
tor /c. This relation applies for the case of bremsstrahlung emis- 
sion, where the cooling function varies with temperature, T, such 
that A oc T^'^^, wh ich renders Ap a function of entropy only 
jPonahue etal]|2005h . 

In the limit of conductive thermal balance, where r = Ap, 
Equation |4] can be rearranged to yield an expression for the corre- 
sponding implied suppression factor: 



/e ~ 62.5 



S3 



(5) 



in terms of the radius, r (in kpc) and the entropy, S (in keV cm ). 
The variation of this quantity with scaled radius is plotted for each 
cluster in Fig.[T3] smoothed using the LOESS function in R and split 
by cool core status. Also shown is the threshold of thermal stabil- 
ity, corresponding to conduction at the Spitzer rate, when fc — 1. 
The striking aspect of Fig. [13] is the clear difference between the 
profiles of CC and non-CC clusters, in terms of both their shapes 
and locations on the plot: CC clusters having higher normalization 
and largely negative gradients compared to non-CC profiles, which 
all lie well within the conductively stable region with mostly pos- 
itive gradients. Once again, the outlier cluster Abell 2142 (as seen 
in Fig. Is} stands out as an anomaly (the only non-CC curve aligned 
from bottom right to top left), which nevertheless retains similar- 
ities to both categories in that it lies within the range of the non- 
CC profiles while bearing closer resemblance to the CC profiles in 
shape. 

The implication of Fig. [13] is that non-CC clusters are cer- 
tainly capable of being stablilized by conduction alone, operating 
at no more than ~10 per cent of the Spitzer rate and even much 
less effectively than that within their inner regions. By contrast, the 
CC clusters require val ues of fc > 0.1 within the peak tempera- 
ture radius (~ 0.15rsnn: |PaDer tt . rising sharply towards the centre 
and in some cases crossing into the region of thermal instability. 
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corresponding to /c > 1 (for NGC 5044, Abell 262, Abell 478, 
Abell 496 and 2A 0335+096). For these five clusters, an additional 
heat source would be required in order to maintain thermal balance. 
With the exception of A478, these are 4 of the 5 coolest CC clus- 
ters (see Table|2}, which is unsurprising, since thermal conduction 
operates much less effectively at lower temperatures as the Spitzer 
conductivity, k oc y^/^. Furthermore, all except Abell 496 show 
evid ence of AGN-relat ed disturbance in the form of X-ray cavi- 
ties teirzan et al1l2004l) and both X-ray and radio disturbance in 
the case of NGC 5044 (David et al, in prep.). These results are in 
agreement with the findings o f lVoigt&Fabia3j2004 . who studied 
four of the same clusters: for Abell 2029 & Abell 1795 they also 
concluded that conductivity at or below the Spitzer level was able 
to balance cooling everywhere inside the cooling radius, whereas 
this was not the case for Abell 478 and 2A 0335+096. 

Theoretical considerations indicate that magnetohydrody- 
namic turbulence in the I CM could give rise t o conduction sup- 
pression factors of ~0.2 dNaravan & Medvedevll200lh . which is 
al so consistent with results from the dire ct numerical simulations 
of iMaron. Chandran & Blackmanl ( |2004|) . who favour fc ~ 0.1- 
0.2. This would imply conductive stability for the non-CC clusters 
at all radii and for all CC clusters outside ~0.05-0.1r5oo (Fig.ll3b. 
Taking all the fc values for the non-CC clusters, we find the spread 
of values to be well fitted by lognormal distribution with a mean log 
(base 10) of -1.50±0.03 and log standard deviation of 0.39±0.02 
corresponding to 0.032 with a la range of 0.013-0.077. 

Why, though, should a threshold for conduction stability give 
rise to a well-separated bimodal distribution of cluste rs, with a 
'zone of avoidance' in between? iDonahue et al. I ( I2OO5I) point out 
that clusters exceeding the threshold for conductive stability will 
continue to cool until some other form of heating intervenes, which 
would therefore lead to diverging populations. In which case, the 
fact that the non-CC fc profiles are themselves not flat suggests 
that this threshold is spatially varying. Since the likely dominant 
cause of conduction suppression is magnetic fields, it follows that 
the configuration of magnetic field lines varies significatly with ra- 
dius. Consequently, it is possible that the variation of fc profiles 
amongst the non-CC clusters reflects the intrinsic variation in in- 
tracluster magnetic field configurations in the cluster population. 
In any case, the non-CC fc profiles (except for A2142) all drop 
sharply inside ~0. Irsoo, reaching values within the typical scale of 
the central cluster galaxy (~0.02r5oo) consistent with the higher 
suppression factors of > O.OI present in the interstellar medium 
of cluster galaxy coronae dSun et alh OOi). This indicates that the 
transition within the ICM to higher suppression factors on galaxy 
scales may occur smoothly, suggesting a gradual variation in mag- 
netic fields. 

By contrast, the presence of cold fronts in the ICM implies 
sharply discontinuous magnetic fields, especially considering that 
such features occur at larger cluster radii, where fc is greater. For 
example, in the case of Abell 2142, detailed analysis of the temper- 
ature gradient across the cold front indicat es that co nduction must 
be suppressed by at a factor of 250-2500 fettori & Fabian 2000), 
which is certainly far below the levels implied in Fig.[T3] although 
these amount to globally-averaged estimates. Such discontinuous 
configurations are possible in this case, since the magnetic field is 
likely to be stretched by tangental gas motions near the cold front, 
making it st ronger and changing its structure comp ared to the rest 
of the ICM jVikhUnin, Markevitch & MurravlbOOlV 

A notable feature of Fig. [T3] is the absence of any sub- 
stantially flat profiles, in between the strongly negative gradients 
of the CC clusters and the mostly strongly positive (except for 



A2142) gradients of the non-CC clusters. A flat profile, lying 
along the conduction-stabilized threshold with a constant value of 

i 'e ^ 1, w ould imply an entropy profile of the form S oc r^^^ 
Donahue et al. 2005), which is close to the cross-over point of 0.69 
between the two distributions of entropy profile logartithic slopes 
found in Section |5] (also plotted as the dashed horizontal line in 
Fig.O. Thus the divergence between CC and non-CC profiles seen 
in Fig.[T3]mirrors the bimodality seen in the slopes of the entropy 
profiles described above. 



6 DISCUSSION 

The above results clearly indicate that bimodality is present in the 
cluster population, structured along the dichotomy between cool 
core and non-cool core clusters. Furthermore, the implication is 
that thermal conduction alone is a plausible mechanism for stabiliz- 
ing the intracluster medium in non-CC clusters, whereas additional 
heating from galaxy feedback is necessary to achieve the sam e 
in CC clusters- a conclusion also reached by IGuo et al] ( l2008h . 
Conductive heat transfer could certainly explain why no signifi- 
cant temperature decline is observed in non- CC clust ers, despite 
the short cooling times of gas in their cores jPaper J) . However, 
by contributing to the heating of the cluster core, conduction also 
acts to reduce the amount of energy input required from additional 
sources such as AGN, in order to maintain thermal stability. This 
may therefore account fo r the observed lack of eviden ce for strong 
AGN heating in clusters jMcNamara & Nulsenll2007l) . 

Nevertheless, while feedback heating may not be necessary to 
stabilize non-CC clusters, it is clear from the results of the entropy- 
temperature scaling analysis in Section 14.2.21 that both CC and 
non-CC clusters show significant departures from self-similarity. 
Such similarity breaking is the unmistakable signature of non- 
gravitational physics, which suggests that the non-CC clusters must 
also have been impacted significantly by feedback (if not also ra- 
diative cooling) in their lif etimes. This possibility is consistent with 
the theoretical model of iMcCarthv et alj 120081) , where non-CC 
clusters are formed from material that has experienced higher lev- 
els of preheating. In this picture, the influence of conduction could 
help to segregate the cluster population and stave off the formation 
of cool cores in the most strongly preheated systems. 

A related aspect is the role of mergers in cluster evolution, par- 
ticularly given the close association between signatures of recent 
disruption, such as radio halos, and the absence of a cool core (e.g. 
lMiIIion& Aiienli2008) . Although suggesting a key role for ther- 
mal conduction in sustaining non-CC clusters, our results neverthe- 
less certainly do not rule out merging as an alternative explanation 
for their existence. While recent simulations have concluded that 
cluste r mergers cannot permanently erase cool cores l lPoole et al.l 
l2006h . it is possible that the additional influence of conduction 
could achieve this outcome and thereby provide an alternative path 
to a non-CC state. In this situation, the temporary fragmentation of 
a cool core that can take place in a merger (Poole et al. 2008) could 
plausibly lead to conductive dissipation of the resulting blobs of 
cool gas, whose size can easily fall below the Field length, pro- 
vided the magnetic field configuration is favourable. Since conduc- 
tion lowers the threshold necessary to transition from a state where 
additional feedback is required (i.e. in a cool core), to one where 
conduction alone can maintain stability (i.e. a non-cool core), the 
transformation from CC to non-CC status is correspondingly more 
achievable. 

However, notwithstanding this appealing explanation of bi- 
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modal populations, it is also clear that the key role that conduc- 
tion can play in stabilizing clusters against cooling instabilities 
cannot easily extend to galaxy groups, where the lower gas tem- 
peratures render it much less effective (the conduction coefficient, 
K oc T^^^). Therefore, the impact of merging activity may be more 
short lived in galaxy groups, without the contribution of significant 
conductive heat transfer to impede re-formation of a cool core. This 
suggests that non-cool core groups are likely to have been recently 
disrupted; by contrast, conduction in clusters could sustain a non- 
CC state in post mergers long after disruption. Nevertheless, it re- 
mains to be seen whether a consistent thermodynamic picture of 
ICM evolution in both groups and clusters can be developed. 

The similarity of the ICM metallicity in CC and non-CC clus- 
ters is noteworthy as the only property of the gas that does not dif- 
fer significantly between the two types, at least within the region of 
overlap seen in this sample. This fact alone provides evidence that 
strong mixing of gas cannot preferentially have affected non-CC 
clusters- as a result of merging activity, for example. However, the 
width of the metallicity peak in the cluster cores appears signifi- 
cantly (> 2x) broader than the stellar distri bution of the brightes t 
cluster galaxy (BCG), as previously noted bv Rebusc o et alj j2005l) . 
which implies some mechanism to transport metals into the ICM. 
iRebusco et alj speculate that this might be caused by AGN-driven 
gas motion, as also favoured by the recent theoretical modelling of 
iRasera et alJfeOOsh . However, if AGN are responsible for diffusing 
the enriched gas then the similarity of cluster metallicity profiles 
implies an equally prominent role for AGN outflows in non-CC 
clusters, despite the lack of evidence for significant AGN disrup- 
tion in such cases. 



7 CONCLUSIONS 

Using th e statistically-sel ected sample of 20 galaxy clusters pre- 
sented inl Sanderson et al . (2006), we have studied the density, en- 
tropy and metallicity of the intracluster medium as a function of 
radius, focussing on the comparison between clusters with cool- 
cores (CC) and those without. We describe an improved method 
of estimating the cluster mean temperature and fiducial scaling ra- 
dius rsoo , which we use to explore systematic trends in cluster gas 
properties across the sample. 

We find that the gas density is systematically higher in the 
cores of CC clusters, and that the ICM is progressively depleted 
in less massive systems. We also find a clear departure from self- 
similar scaling in the gas entropy which is consistent with the modi- 
fied scaling of 5" ex T*" with 6 = 0.65 from lPonman et al.l ( l2003l) : at 
O.lrsoo the best-fit scaling is 6 = 0.66 ± 0.10 and 0.71 ± 0.21 for 
CC and non-CC clusters, respectively. However, the dependency 
on temperature strenghens when all the clusters are combined, to 
a nearly self similar value of fe = 0.92 ± 0.12, and similar results 
are obtained in all three cases for entropy measured at O.lSrsoo 
versus mean temperature. This demonstrates that similarity break- 
ing (i.e. S (X y~2/3-j exist in the separate populations of CC 
and non-CC clusters, even while the combined population shows 
consistency with self-similarity. 

The metallicity of the gas shows no evidence of a system- 
atic variation with T, but declines with radius such that Z cx 
^-o.3i±o.04 Qyj^j^jg o.014r5oo (comparable to the size of any cen- 
tral dominant galaxy), for both CC and non-CC clusters alike. In- 
side this point there is substantial divergence in the metallicity, with 
a few CC clusters showing sharply decreasing Z towards the cen- 
tre while others possess continually rising profiles. At large radii. 



there is no indication of any flattening in the metallicity profile to 
at least ~0.5r5oo, where the Chandra field-of-view limits the data. 
We study gas metallicity as a function of entropy and find a strik- 
ing lack of low-metallicity gas (i.e. < 0.4-0.5 Solar) with low en- 
tropy {S < 200keVcm^). Above ~100keVcm^ the metallicity 
declines with increasing entropy in an identical fashion for both 
CC and non-CC clusters. 

We address the issue of bimodality in cluster properties by 
studing the distribution of logarithmic slopes obtained from power- 
law fits to individual cluster entropy profiles (i.e. S oc r"). We find 
that a double-Gaussian distribution is strongly preferred over a uni- 
modal Gaussian distribution, using maximum likelihood fits to the 
unbinned values, employing both the Bayesian and Akaike Infor- 
mation Criteria model selection tests. The best fitting means of the 
two distributions are a = 0.92 ± 0.04 and 0.48 ± 0.04, with a stan- 
dard deviations of 0.1 in both cases. Given the statistically selected 
nature of the sample, this demonstrates that two distinct categories 
of cluster exist, which has important implications for models of 
galaxy feedback and cluster similarity breaking. 

We explore the impact of thermal conduction on the ICM by 
studying the implied conduction suppression factor, /c, as a func- 
tion of radius. We find that the profiles of fc differ sharply between 
CC and non-CC clusters consistent with two distinct populations. 
We conclude that conduction alone is capable of stabilizing non- 
CC clusters against catastrophic cooling, while in CC clusters some 
feedback is required in addition to conduct ion to main tain thermal 
balance, in agreement with the findings of iGuo et alj (ZOOS). Tak- 
ing all the fc values for the non-CC clusters, we find the spread of 
values to be well fitted by lognormal distribution with a mean log 
(base 10) of -1.50±0.03 and log standard deviation of 0.39±0.02 
corresponding to 0.032 with a la range of 0.013-0.077. 
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